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1. Introduction 

To tile means to cover a given part of space, without overlaps, using a small number, 
often just one, of different types of objects. When we are covering the floor of a room 
using identical rectangular tiles (in the ordinary sense of the word) we lay the copies of 
the tiles in a regular way next to each other, leaving no gaps. 

Your floor can be as boring as the rectangular floor on the left or as interesting as the 
Escher lizard tiling on the right: 




In both these floor tilings there is only one tile used, a square tile on the left and a lizard 
shape tile 
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on the right. 

Beyond aesthetics and decoration of the tile there is another difference between these 
two examples which is very significant from our point of view. To fill the floor using the 
square tile we only need to translate the tile around; we never have to turn it. This is not 
the case with the lizard tile. To fill the floor the various copies of the lizard tiles have to 
be turned appropriately so as to fit perfectly with each other. 

It is with the seemingly boring case that we shall occupy ourselves in this survey: the 
single tile that we have at our disposal is only allowed to translate in space. We cannot 
rotate it or reflect it. True, in the vast majority of examples in the literature (see e.g. 
[13]) that exhibit interesting behavior (such as undecidability or aperiodicity, on which 
we'll have more to say later) the allowed group of motions is usually the full group of rigid 
motions in space and therefore rotations are allowed (and some times even reflections as 
well) . One would indeed be hard pressed to find eye-catching examples such as the Penrose 
"kite and dart" aperiodic set of tiles shown below. 




But we have good reasons for restricting ourselves to tiling by translation which we hope 
to expose in this survey. To begin with let us state that tiling by translation is still full of 
unresolved questions with connections to number theory, Fourier analysis and the theory 
of computation, and that happens even in dimension 1. If it is hard to imagine how one 
can seriously study questions of tiling in dimension 1 let us point out that the tiles need 
not have a "nice" shape. For instance, one can easily show by trial and error that the set 
E = {0, 2, 3, 5} (show below in red) 

• o • • o • oooooo 

can tile the integers Z by translation only and that such a tiling has period 12. This trial 
and error process can become slower and slower if one considers larger finite sets of integers 
and simply asks if they could be used to tile Z. In fact, it is not known how to do this in 
essentialy any other way. This reflects how little we understand this tiling phenomenon. 
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Here is another example, of a more geometric flavor, to convince the reader that taking 
away the freedom to rotate and reflect the tile does not take away the fun. For this consider 
the so called notched cube (left) 




and try to visualize a way to fill 3 dimensional space with it (here it would not help to rotate 
or reflect it). We are making no assumptions here about the lengths of the rectangular cut 
(notch) at one corner of the cube or rectangle. Whatever these lengths are this notched 
cube can indeed tile space by translation [18]. Yet the simplest way we know how to prove 
this is using Fourier analysis, making absolutely no use of geometric intuition. Same holds 
for the extended cube (above right) and this seems to be even harder to visualize. At least 
the notched cube "forces" a corner of a copy of the tile to fill the notch in another copy 
and this can at least get you started but no such restriction is obvious for the extended 
cube. 

By now we should fix the rules of the game. What does it mean for a subset Q of W 1 to 
tile when translated at the locations A C 18L d (another discrete set)? Simply that the 
copies 

Q + A, A 6 A, 

are mutually non- overlapping, and cover the whole space. Here, as is customary, we denote 
by O + A the set 

{lu + A : we^}, 




or, in other words, the translate of f2 by the vector A. And what does it mean for two 
sets to be non-overlapping? One could, for instance, demand that the interiors, but not 
the boundaries, of the two sets are disjoint. But it turns out that the most convenient way 
is to demand that the intersection of these two sets has zero volume (or area in dimension 
d = 2, or length in dimension d = 1) or, more precisely, zero Lebesgue measure. 

2. Using the Fourier Transform 

2.1. Tiling in the Fourier domain. Having specified what we mean we turn now to 
exploiting the technology that we know, and in this case we find it very fruitful to define 
tiling simply by the equation 

(1) Xfi( x ~ A) = 1, (almost everywhere for x € 

xeA 
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We now have an equation! as a physicist would exclaim. A mathematician should never 
underestimate the power of formal manipulation and the phenomena it could reveal. So, 
let us rewrite our equation (ignoring from now the exception of measure): 

(2) xn * <5a = 1- 

Here the * operator denotes convolution and 5\ is the measure 




<5a = ^ 8\, 

AsA 

where 5\ is a unit point mass sitting at point A € M. d . (We recall that / * 5\(x) — f(x — A), 
so that convolving a function / with 5\ merely translates the function by the vector A.) 
The object 5\ is a compact way to encode all information that's contained in the set of 
translates A in a way that we can operate algebraically on it. 

An equation such as (2) of course begs to be subjected to the Fourier Transform 1 as 
the Fourier Transform (FT) behaves so nicely with convolution (the FT of a convolution 
is the pointwise product of the FTs of the convolution factors: fg = fg). The definition 
of the FT / of a function / that we use is 

(3) /(£)=/ f(x)e- 2 ^*dx, (eG< / |/| < oo). 

Jm. d Jm. d 

This definition is sufficient for the function / = xn which is integrable but one needs to 
define the FT of a measure by duality (see for instance [22] for the details). Postponing 
such worries till later we take the FT of both sides of (2) to get 

(4) • 5 A = S . 

The measure So on the right is just a point mass at point 0. The support (where it is 
"nonzero") of this is {0} and therefore that must be the support of the left hand side of 

(4) . Since that is a product it follows that, apart from point 0, wherever xn is nonzero, S\ 
must be zero so as to kill the product. Let us summarize this in the following necessary 
(and often sufficient) condition for tiling: 

(5) supp^C {xh = 0}U{0}. 

A formal proof of (5) can be looked up in [22]. Let us make no comment here about when 
this condition is sufficient. But we must point out that the object S\ is not necessarily a 
measure but a so called tempered distribution [35]. 

2.2. A problem from the Scottish Book. Nice reformulation but what could this buy 
for us? 

Following the time-honored tradition of our trade, let us first generalize: whatever we've 
said so far holds just as well for an (almost) arbitrary function / (nonnegative, integrable 
would do) in place of xn- In other words, if one has 

(6) ^ f(x - X) = £ = Const., 

AeA 

for almost every x 6 K d then it follows that 

(7) supp^C Z(/)U{0} where Z(/) = {/ = o}. 



This of course depends on what training you have been subjected to. 
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Whenever (6) holds we will say from now on that / tiles space with A at level £. And 
Steinhaus suggested that perhaps f(x) = e~ x has this property (although he quickly 
disproved it). 

Let us apply our new tool to Problem No 181 in the famous Scottish Book 2 [31], a 
problem posed by H. Steinhaus and partly solved by him around 1939. 



(in the sense of H. Steinhaus) for every couple f,,/,(f, * U)l 


181 


Find a continuous function {or perhaps an analytic one) 
f{x), positive and such that one has 

OP 


I V ■ 

H. STEINHAUS 


E f{X + B) ■ 1 
ft - - » 




(identically in x in the interval -a> < x < + »); examine 
whether {\/-^e~'' is such a function; or else prove the 
impossibility; or else prove uniqueness. 




Addendum. The function (l/VVJe"*' does not have the 
property — this follows from the sign of the second derivative 
For v = of the expression 




E -U-'-"' 1 . 

- v H. Steinhaus 





In our language the question essentially was if there is a positive, continuous and even 
analytic function / that tiles M with Z. 

2 

Whoever has been exposed even a little to the FT knows that e~ x has the nice property 
that its FT transform is essentialy itself (times a constant and possibly rescaled, depending 
on the normalizations used in the definition of the FT) and therefore it has no zeros at 
all and therefore (7) has no chance of holding as <5a has no place to be supported at apart 
from 0, and it's not a constant. And this argument does not use the fact that A = Z so 

2 

e~ x does not tile M with any set of translates. 

An analytic solution to this question of Steinhaus can also be provided using our tech- 
nology. First of all, let us make the important remark that 

(8) 6 Z d = 6 X d . 

This is the deservedly famous Poisson Summation Formula [36, Ch. 9] 

n n 

in disguise. In this case, when 6a is a measure, condition (7) is also sufficient for tiling [22], 
so it is sufficient to find a positive function / whose FT vanishes on Z \ {0}. Analyticity 
of / will be guaranteed if / has compact support. It is then easy to check that one can 
take as / the sum of two triangle functions with incommensurable base lengths, 



as the FT of each such triangle is a nonnegative function (a Fejer kernel) whose zeros are 
at integer multiples of the reciprocal of its half-base length. Since these will never match 

2 A collection of problems from the Scottish Cafe in Lwow (then in Poland now Lviv in the Ukraine) 
where a number of important mathematicians met for a number of years exchanging problems. These 
problems and whatever solutions were recorded in a notebook that formed the basis for [31]. 
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for the two Fejer kernels their sum is always positive. Since we also take the bases of the 
two triangles to be supported in ( — 1, 1) it follows that (7) holds and we have tiling. 

2.3. Filling a box with two types of bricks. Let us give another amusing application 
of the FT method to tiling problems. This will be a Fourier analytic proof [21] of a result 
of Bower and Michael [5] . Suppose that you have two types of rectangular bricks at your 
disposal, type A of dimensions a\ x (12 and type B of dimensions b\ x 62 (we're stating 
everything in dimension 2 but everything works in any dimension) and your task is to tile 
a rectangular box, say Q — (—1/2, 1/2) 2 , with copies of bricks A and B. The bricks may 
be translated but not rotated. 

We will show that this is possible if and only if you can cut Q along the x or along the 
y direction into two rectangles each of which can be tiled using bricks of one type only. 
A generic tiling of the type shown below left implies the existence of a separated tiling as 
shown below right. 





In other words, having two types of bricks at your disposal does not demand much ingenuity 
on your side in order to exploit for tiling purposes. 

A simple calculation shows that ifC = ( — y, y) x (— %) is a centered c\ x C2 box 



then 



3in(7rci£) sin(7TC2?7) 



and therefore xc vanishes at those points, shown below, 











V 
































l/c 2 
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where the £ or the r/ coordinate is a non-zero multiple of 1/ci or I/C2 respectively 

Suppose now that we can tile Q by translating copies of brick A to the locations T and 
copies of brick B to the locations S. In other words 



(9) 

which we rewrite as 



XQ = XA* St + Xb * Sg 
and take the FT of both sides to get 

(10) xq(£, v) = <h(€, v)xa(£, v) + 0s(£. v)xb(^ v), 
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where (fix = <5r and <f>S — °S are two trigonometric polynomials. Since Q — (—1/2, 1/2) 2 
we have 

(11) z(xq) = {x^ = o} = {Zez\{o}or V ez\ {o}}. 

But, because of (10), \Q must vanish on the common zeros of \X and xb, for instance at 
the points (1/a.i, I/&2) and (l/&i, 1/0-2)) which implies, because of (11) 

(12) (l/oi 6 Z or l/6 2 6 Z) and (l/&i G Z or l/o 2 G Z). 

If (12) is true because l/oi, l/a 2 6 Z then brick A alone can fill Q. Similarly if (12) is 
satisfied with l/&i, I/&2 6 Z then brick B alone suffices. 

What happens if l/oi, 6 Z? Since we have assumed a tiling of Q in (9) it follows, 
by traversing the y-axis, that there are nonnegative integers k, I such that 

1 = ka 2 + lb 2 . 

Cut now the Q box parallel to the x-axis at height ka 2 from the bottom as show here: 



1 













ka 2 



Now it is clear that brick A can fill the lower box (since 1 ja\ G Z) and brick B can fill the 
upper box (since l/&i G Z). The remaining case l/a.2, I/&2 G Z is treated similarly. 

3. Discrete tilings 

3.1. Tilings of the integers and periodicity. Let us now focus on tiling the integers. 
Let A C Z be a finite set and A G Z. We say that A tiles Z with A at level I if the copies 
A + A, A G A cover every integer exactly I times. In other words 

xaOe — A) — £, for all x G Z. 

AeA 

We denote this situation as A + A — PL. 

We say that a tiling is periodic with period t G Z if A + 1 — A. It is a basic fact proved 
by Newman [32] that all tilings of the integers at level 1 are periodic. 

Indeed, suppose that A = {0 — a\ < . . . < a^} is a finite set (we may freely translate A 
without changing the tiling or the periodicity property so we assume it starts at 0) and 
A + A — Z is a tiling at level 1. Fix any x G Z and write W x = x, x + 1, . . . , x + — 1 
for the "window" of width afc — 1 (one less than A) starting at x. We claim that the set 
A fl W x determines A. Let us show that it determines A to the right of x + — 1. It will 
determine A to the left of x by the same argument. 

It is enough to decide, looking at A n W only, if x + cik G A or not. (We then repeat for 
x + dk + 1 and so on.) Observe that for any A G A fl (—00, x) the set A + A is contained 
in (—00, x + dk — 1], so any such copy cannot be used to cover x + a*;. Clearly no copy of 
the form A + A for A > x + at can be used for that purpose too. We conclude that x + 
is covered by some copy A + A with A G W U {x + a^}. Inspecting A fl W x we can tell if 
the relevant A is in W x or not. If it is in W x then x + £ A since this would lead to the 
copies A + A and A + (x + at) to overlap at x + a^. If it is not in W x then necessarily 
x + afc G A, and this concludes the proof of the claim. 

How many different values can the set A fl W x take? Clearly it can take at most 2 afc 
different values as there are two choices (in A or not in A) for each x G W x . This means 
that there are two different x, y G {0, 1, . . . , 2 afc } for which A(~)W X is a translate of AnW y . 
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It follows that A + (y — x) = A. We have proved that every tiling has a period which is at 
most 2 D where D is the diameter of the tile. 

There is a similar result for tilings of the continuous line by translates of a function 
[23, 28, 22]. Combinatorial arguments do not seem to be enough here and the Fourier 
analytic technology along with some deep results of Harmonic Analysis are being used for 
the proof. 



3.2. Tilings of the finite cyclic group. The fact that a tiling of the integers is periodic 
allows us to view it as a tiling on a smaller structure, a cyclic group. Indeed, assume that 
the tiling A + A = Z has period n, that is A + n — A. Define the set 

A — A mod n C {0, 1, . . . , n — 1} 

by taking for each A £ A its residue mod n. It follows from the n-periodicity of A that 
A — A + nZ (this is again a tiling or a direct sum: every element of A can be written in a 
unique way as an element of A plus an element of nZ) . Hence we have 

Z = A + K = A + K + nL 

with all sums being direct. Taking quotients we obtain that the cyclic group Z„ — Z/(nZ) 
can be written as a direct sum (tiling) 

(13) Z„ = A + A. 



In this case we obviously have n — \A\ 



A 



Let us stop here to make two side remarks about periodicity. The first is that tilings 
of the cyclic group Z„ can also be periodic. Indeed, assume Z n — A + A and there exists 
/ k £ Z„ such that A is periodic by k, i.e. A + k — A. Then we can reduce the set 
A modulo k, A = A mod k C {0, 1, . . . , k — 1}, and conclude as above that Z& — A + A. 
Therefore, periodic tilings of Z n can be regarded as tilings of a smaller group (it is 
trivial to see that k automatically divides n). It is thus natural to ask whether certain 
cyclic groups Z n admit only periodic tilings, i.e. whenever Z n = A + A, then either A or A 
must be periodic. These were called "good" groups by Hajos [14] (the notion also makes 
sense in the more general setting of finite Abelian groups, not only cyclic groups). It turns 
out that some groups indeed have this property, and Sands completed the classification 
of good groups in [38, 39]. In particular, the good groups that are cyclic are Z„ where n 
divides one of pqrs, p 2 qr, p 2 q 2 or p m q, where p, q, r, s are any distinct primes. The cyclic 
group of the smallest order which is not good is Z72. 

The other remark concerns the connection of tilings of cyclic groups to music com- 
position. As explained above, periodic tilings of Z„ are mathematically speaking less 
interesting because they can be considered as tilings of some smaller group Z&. It turns 
out that non-periodic tilings of a cyclic group Z„ are also more interesting from an aes- 
thetic point of view, and they are called Vuza- canons in the musical community. The 
interaction of mathematical theory and musical background has been extensively studied 
in recent years [3, 1, 2, 10, 42]. Finding all non-periodic tilings of Z„ is therefore motivated 
by contemporary music compositions. Of course, the problem makes sense only if Z™ is 
not a good group (otherwise all tilings are periodic). Fripertinger [9] achieved this task for 
n = 72, 108 while the authors [26] gave an efficient algorithm to settle the case n — 144 
(the algorithm is likely to work for other values like n — 120, 180, 200, 216, beyond which 
the task simply seems hopeless.) 

Let us now return to the analysis of tilings of the cyclic group Z„. The Fourier condition 
(5) for tiling takes exactly the same form here and is much simpler to prove as no subtle 
analysis is required (no integrals, only finite sums are involved). The Fourier transform of 
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a function / : Z n — > C is denned as the function / : Z n — > C given by 

n-l 
3=0 

where £„ — e 2m / n is a primitive n-th root of unity. 

In the finite case the roles of the tile A and the set of translations A in (13) can be 
interchanged, so let us adopt a more symmetric notation 

Z n = A + B, 

in which A, B C Z„ are just two subsets, necessarily satisfying \ A\ ■ \B\ — n, and such that 
every element of Z„ can be written uniquely as a sum of an element of A and en element 
of B. The Fourier condition now takes the form 

(14) Z n = {0}UZ(x2)UZ(x^). 

As we saw in the section §3.1 every tiling of Z has a period which is at most 2 D , where D 
is the diameter of the tile. We will see in this section how the easy combinatorial argument 
of §3.1 can be replaced with an argument that has the Fourier condition (14) as a starting 
point, uses some well-known number-theoretic facts (cyclotomic polynomials) and gives 
much better results. We will describe the result in [20] and by I. Ruzsa in an appendix in 
[41] since it is simpler than the current best results in [4]. 

So, suppose that A C {0, 1, ... , D} is a set of integers of diameter < D and that A tiles 
Z with period M. This implies that A (to be precise, A reduced mod M) tiles the cyclic 
group Z M 

Z M = A + B. 

Assume also that M is the least period. This implies that if g £ %u is such that B = B + g 
then g = 0. (Otherwise the original tiling would be periodic with period g.) 

Now we use the Fourier transform without really mentioning it, using instead polynomial 
terminology. Let 

A(x) = J2* a , B(x) = Y,x\ 

aeA beB 

be the polynomials defined by A and B. It is very easy to see that the tiling condition 
Z M — A + B is the same as 

A{x)B{x) = 1 + x + x 2 + • • • x M_1 mod x M - 1, 

or, in other words, that 

r M _ i 

(15) x M - 1 divides A(x)B(x) , 

x — 1 

and this implies that all M-th roots of unity except 1 are roots of the product A{x)B(x). 
(This last statement is the equivalent of (14) in algebraic language.) 

The M-th roots of unity ( 3 M = e 2m i/M ^ = 0, 1, . . . , M — 1, are grouped into cyclotomic 
classes. Two roots C,° M , C,\j belong to the same cyclotomic class if and only if the greatest 
common divisors (j, M) and (k, M) are the same. If a root of unity in a certain cyclotomic 
class is a root of a polynomial with integer or rational coefficients then so are all the other 
roots in the same cyclotomic class. If d is a positive integer then the <i-th cyclotomic 
polynomial is the monic polynomial which has as roots all the primitive d-th. roots of unity 

*d[x)= J] (x-e 2 ^ d ). 
l<j<d,(j,d)=l 
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We have deg^^x) — (f>(d) (the Euler function) and we will use the estimate [15, p. 267] 

(16) C r^~^ ^ ^ ^ d ' 

log log a 

which holds for some constant C and all large enough d. It turns out that each is an 
inrreducible polynomial with integer coefficients and one can write 

x m - 1 = n 

d divides M 

Coupled with (15) this implies that if d > 1 is a divisor of M then ^(x) divides the 
product A{x)B(x). 

Let now & si (x), . . . , <& Sk (x) be all cyclotomic polynomials <3> s (:r) with s > 1 that divide 
A(x), written once each and numbered so that 1 < si < s 2 < ■ ■ ■ < Sfc. Since deg $ s — <f>(s) 
it follows that 

(17) 0(si) + • • • + 0(s fc ) < deg < D. 
We have 



X! Si - C^2<j){si) log log s t (from (16)) 



i=l i=l 
fe 

<C^</>( Sl )loglogD 

i=l 

<CL>loglogD (from (17)). 

From this inequality and using the fact that the Sj are different integers it follows that 
(18) k < Cy?D log log D, 

as obviously the worst case is if the Sj are k consecutive integers and in that case their 
sum is > k 2 . 

Define now t = Yli=i s «' so tnat au cyclotomic polynomials that divide A(x) are also 
divisors of x l — 1. It follows that 



,M 



x 



1 divides {x l - l)B(x) 



and this means precisely that B has period t, hence t > M as we assumed M to be 
minimal. Using (18) and the bound Sj = 0(D 2 ), for instance, we get the following bound 
for the least period M: 

(19) M <t< exp(C/DlogD^loglogL»), 

for an appropriate constant C. This is a much better bound, in terms of its dependence 
on D, that the bound M < 2 that we obtained with combinatorics alone. 

3.3. Periodicity in two dimensions and computability. In contrast to the one- 
dimensional case, where every tiling by a finite subset of Z is periodic, it is very easy 
to see that in Z 2 there are tilings which are not periodic. One has to be a little careful 
with periodicity in dimension 2 and higher though. We shall call a set A C 7L d periodic 
if there exists a full lattice of periods, i.e. if there exist d linearly independent vectors 
ui,...,Ud which are all periods of the set 

A + u\ — A + u 2 = ■ ■ ■ = A + u d — A. 

If that happens then every integer linear combination of the Uj is also a period or, in other 
words, the set A has a full-dimensional lattice of periods 

A = span z {u l ,u 2 ,...,u d }. 



TILINGS BY TRANSLATION 



11 



Thus in dimension 2 a set may have a period but not two linearly independent periods, in 
which case we do not call it periodic. One example is the set Z x {0} C Z 2 . A tiling of Z 2 
by a finite subset which is not periodic is very easy to construct. 



On takes A to be a square, for instance the set {(0,0), (1,0), (1,0), (1, 1)} and perturbs 
the usual tiling of it by shifting one column only up by one unit. This destroys the period 
(2, 0) of the tiling along the x-axis, but leaves the period (0, 2) along the j/-axis intact. 

It is not entirely obvious how to construct a tiling of Z 2 which has no period at all, 
but it can be done and let us briefly describe how. The key idea is to start again with 
the usual square tiling and do something to it so as to destroy all periods, not just the 
periods along one axis. The way to do this is to simultanesouly shift a horizontal and 
a vertical column. This is not of course possible with the square tile we used before, as 
simple experimentation will convince you. But it can be done if we "interleave" four tilings 
of the type shown above. Take our tile to be the set 

A = {(0,0), (2,0), (0,2), (2,2)} 

which is what we had before only dilated by a factor of two. Now tile (2Z) 2 using A in 
the usual way, that is by translating A to the locations (4Z) 2 . Using this tiling one can 
get a tiling of Z 2 by tiling the four sets (cosets of the subgroup (2Z) 2 in Z 2 ) 

(2Z) 2 , (2Z) 2 + (1,0), (2Z) 2 + (0,1), (2Z) 2 + (1,1) 

in exactly the same way. But the smart way to do it is to use the nice tiling for the first two 
of the above four cosets, then use, for tiling (2Z) 2 + (0, 1), a tiling like that shown in the 
figure above (which destroys the horizontal periods) and finally use for tiling (2Z) 2 + (1, 1) 
a similar tiling which destroys the vertical periods. 

Having established that in two dimensions there are tilings with no periods at all, let us 
now remark that it is very different to ask for tiles which are aperiodic, i.e. for tiles which 
can tile but only in a way that admits no periods. In fact, the answer to this question is 
not known if we insist that only translations are allowed. 

Conjecture 3.1 (Lagarias and Wang [28]). If a finite subset A C Z 2 can tile by translation 
then it can also tile in a periodic way. 

Let us point out that if more freedom than translation is allowed then tiles (or sets of 
tiles) are known which can only tile aperiodically, so the above conjecture, if true, would 
mean that restricting the allowed motions to translations makes a big difference in this 
respect. 

Let us also remark that it can be proved, at least in dimension 2 [33], that if a set 
admits a tiling with just one period then it also admits a (fully) periodic tiling. 

The property of periodicity is of great importance to the issue of computability, namely 
to whether a computer can decide if a finite set A admits tiling by translation. In still 
more words, we are interested to know if there is a computer algorithm (a Turing machine 
for purists) which, given a finite subset A oil? 1 will decide in finite time if there exists a 
tiling complement of A or not, i.e. if there exists B C Z 2 such that A + B = Z 2 is a tiling. 

We do not require our algorithm to: 

• run fast (only to finish at some point) 
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• find such a tiling complement B if it exists. 

Especially for the second point notice that it would not make sense to want to know B as 
such a set is an infinite object and there is no a priori reason for it to be describable in 
some finite manner. 

It is easy to give an algorithm that would answer NO if A is not a tile but would 
run forever if A is a tile. Here is how this can be done. For each n — 1,2,..., decide 
by trial and error if there is a finite set B such that A + B covers the (discrete) square 
Q n = [— n,n] 2 n Z 2 in a non-overlapping way. There is a finite number of such sets B to 
try as it does not matter what B is far away from Q n . If such a B exists move on to 
the next n. If not then declare A a non-tile and stop. The correctness of this algorithm 
(i.e., the fact that the algorithm will stop for any non-tile A) follows from a fairly simple 
diagonal argument: if A can cover in a nonoverlapping way every Q n then it can also tile 
Z 2 . We invite the reader to think over why this is so. 

The hard thing of course is to have an algorithm that always halts and say YES if A is 
a tile and NO if it is not a tile. Let us point out here that in various other tiling situations 
such an algorithm cannot exist. For instance [12] there cannot exist an algorithm which, 
given a finite collection of subsets of Z 2 , decides if these can tile the plane by translation. 
This remains true even if the number of subsets remains bounded (but large) . 

The connection with periodicity [34] is that the periodic tiling conjecture above implies 
decidability. Indeed, let us assume that any set that tiles can also tile periodically. (We 
emphasize here that we are making no assumption about the size of the periods or their 
dependence on the tile.) This means that it is equivalent to decide if a given set A has 
periodic tilings. 

Let us assume that A C [0, D] 2 . For n — 1,2,..., we can clearly find all subsets B of 
[—n — 2D, n + 2D] 2 such that A + B is a nonoverlapping covering of Q n — [— n, n] 2 . We can 
do this, slowly but surely, by just trying all eligible subsets. If we find none then clearly 
our set is not a tile, we say NO, and we stop. If we do find some we look among them to 
find a set B which "can be extended periodically" . How does such a set look like? 

If an infinite set B C Z 2 is periodic with linearly independent period vectors U\,U2 6 Z 2 
then we can always find two vectors 

ui = (a, 0), u 2 = (0, b), a, b £ {1, 2, 3, . . .}, 

which are also periods (in other words, every lattice in Z 2 contains a sublattice generated 
by two vectors on the x- and y-axes). The mapping (x, y) — > (x mod a, y mod b) maps B 
to a set 

B' C {0, 1, . . . , a - 1} x {0, 1, . . . , b - 1} 
and the periodicity of B means exactly that 

B = (aZ) x (bZ) + B' 

is a direct sum, i.e. every element of B can be written uniquely in the form (am, bn) + b' , 
where m, n 6 Z, b' 6 B' . 

Suppose now that our algorithm looks at a finite set B which is a finite part of a periodic 
set B. This finite part of B arises if we keep those elements of B which are used in a non- 
overlapping covering of Q n (shown as a dashed green rectangle). If n is large enough then 
the set B will look something like this: 



TILINGS BY TRANSLATION 



13 



That is, the set B will contain several "full periods" (the dashed blue a x b rectangles) 
B plus a few incomplete periods near the border of Q n . Suppose now that we have a k x k 
block of full periods in the set B and that k is large enough that both ka and kb are much 
larger than the diameter of A, say larger than 10D. Then we can just extend this block 
periodically (with periods (a, 0) and (0, b)) and we will get a periodic tiling complement 
of A. The reason is that (a) no overlaps arise in this manner or they would have shown 
up in the k x k block already and that (b) we have covering of everything as clearly the 
k x k block suffices to cover the square of side 5D with the same center. 



4. Tiles, spectral sets and complex Hadamard matrices 

In this last section we introduce the notion of spectral sets and describe how it is related 
to tilings and complex Hadamard matrices. For the sake of simplicity, we will restrict our 
attention to finite groups, in particular to cyclic groups Z ra and its powers Z^. 

Although we have seen in the sections above that many fundamental problems about 
translational tilings remain open, the reader will agree with us that the notion of tiling 
is very intuitive and easy to grasp. To the contrary, the definition of spectral sets is 
somewhat more abstract. In order to make it more down-to-earth, we will not use the 
standard definition here, but first introduce complex Hadamard matrices and then use 
them to define spectral sets. 

What are Hadamard matrices? Classically, they are square matrices consisting of el- 
ements ±1 only, such that the rows (and thus also the columns) are orthogonal to each 
other. Complex Hadamard matrices are a natural generalization of this concept. A k x k 
matrix H is a complex Hadamard matrix if its entries are complex numbers of modulus 
1, and the rows (and thus the columns) are orthogonal. Recall here that orthogonality is 
understood with respect to the complex scalar product, i.e. you need to conjugate in one 
of the components, the scalar product of (z\, . . . Zk) and (m, . . . Uk) being Ylj=i ~Zj u j- 

It is well-known that a k x k real Hadamard matrix can only exist if k is divisible by 4. 
On the contrary, complex Hadamard matrices exist in all dimension. Indeed the Fourier 
matrix Fk defined as 

11 1 

1 Ck 



F k := 



/■k-1 
^k 
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where (jfc — e 2m / k , is a complex Hadamard matrix for every k. 

Each element hj_ m of a k x k complex Hadamard matrix H is of the form e 2 ™ p ^ m where 
Pj,m € [0, 1)- We will call the real k x k matrix R formed by the angles p^ m the logarithm 
of H, in notation, R = \og{H) or H = exp(i?). For a real matrix R such that exp(i?) is 
complex Hadamard, we will use the terminology that R is a log- Hadamard matrix. 

Complex Hadamard matrices play an important role in quantum information theory, 
in particular in the construction of teleportation and dense coding schemes [43] , among 
other things. An online catalogue of all known families of complex Hadamard matrices is 
available at [6]. 

Let us now turn to the definition of spectral sets. We identify elements of Z^ with column 
vectors, each coordinate being in the range {0, 1, . . . , n — 1}. Thus, a set S C with r 
elements can be identified with andxr matrix, the columns of which are the elements of 
S (the order of elements does not matter) . We will abuse notation and denote this matrix 

also by S. We also identify Z^ with row vectors whose coordinates are in {0, 1, . . . , n — 1}. 
Accordingly, we will use the notation Z^ — Z^ T (the T meaning transposition). Two sets 
S C Z^ and Q C Z^ T , both of them having r elements, are called a spectral pair if the 
matrix ^QS is an r x r log-Hadamard matrix. In this case Q is called the spectrum of S, 
and S is called a spectral set. We remark that S T Q T = (QS) T (where S T and Q T denote 
the transposed matrices) is then automatically also a log-Hadamard matrix, so we obtain 
that S is a spectrum of Q, which justifies the symmetry of the terminology "spectral pair" . 

For example, the whole group Z n = {0, 1, . . . , n — 1} is spectral, and its spectrum is the 
whole dual group Z^, giving rise to the logarithm of the Fourier matrix, log(F n ) = -Z^Z n . 

We saw that the tiling Z„ — A + B is characterized by \A\ ■ \B\ = n and the Fourier 
condition Z n = {0} U Z (xa) U Z (xb) in equation (14). After some inspection one sees that 
spectral sets also admit a characterization by Fourier analysis. Indeed, using the definition 
above we obtain that S C Z^, \S\ = r, is spectral if and only if there exists Q C Z^ T , such 
that |Q| = r and 

(20) qj-q k eZ(x~s) 

for each qj / 6 Q- (This is an equivalent way of saying that the jth and kth rows of 
the matrix exp(iQS') are orthogonal.) 

What are the connections of tiles and spectral sets? Notice that by the orthogonality 
conditions the matrix U — ejcp(^QS) is unitary. This implies the following important 

fact (for a formal proof, see [25]): 

(21) ^|xS| 2 (*- g ) = |S| 2 = r 2 , 

q&Q 

for all x G Z^ T . In other words, the function \xs\ 2 tiles Z^ at level \S\ 2 when translated 

at the locations Q. One can also say that the set Q tiles the group Z^ with the weighted 
translates defined by the nonnegative function 




The question thus arises naturally: does Q tile Z^ in the ordinary sense, i.e. does there 

exist a set P c Z^ such that Q + P = Z^? A famous conjecture of Fuglede [11] concerns 
exactly this scenario: it states that a set S is spectral if and only if it tiles the group 
(the conjecture was originally stated in the Euclidean space R d but it makes sense in any 
abelian group). 

Several positive partial results were proven concerning special cases of Fuglede's conjec- 
ture (see e.g. [16, 17, 19, 27]), before T. Tao [40] showed an example of a spectral set in 
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M 5 which did not tile the space. Tao's example was based on considerations in the finite 
group Zjj. It also implies the existence of such sets in any dimension d > 5. Subsequently, 
counterexamples in lower dimensions (4 and 3, respectively) were found by the authors 
[29, 24]. All these examples are based on considerations in finite groups Z^ and, ultimately, 
on the existence of certain complex Hadamard matrices. 

The other direction of Fuglede's conjecture could not be settled by Tao's arguments. 
Finally a non-spectral tile in IR 5 was exhibited by the authors [25] by a tricky duality argu- 
ment, also based on considerations in finite groups. Later, counterexamples in dimensions 
4 and 3, respectively, were found in [8, 7]. As of today, Fuglede's conjecture is still open 
in both directions in dimensions 1 and 2. 

As a final remark let us mention another interesting connection of tilings and complex 
Hadamard matrices. As discussed above, Fuglede's conjecture is not true in general, but it 
is true in many particular cases. What could this buy for us? We saw from the definition 
that spectral sets are directly related to complex Hadamard matrices, and the latter are 
very useful in another branch of mathematics: quantum information theory. Could we use 
the connection of tilings to spectral sets (which exists in many special cases) to construct 
new families of complex Hadamard matrices? It turns out that the answer is positive. It 
turns out [30] that "trivial" tiling constructions unfortunately lead to well-known families 
of complex Hadamard matrices (so-called Dita- families) . However, a non-standard tiling 
construction of Szabo [37] was used in [30] to produce previously unknown families of 
complex Hadamard matrices in dimensions 8, 12 and 16. 
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1. Introduction 

To tile means to cover a given part of space, without overlaps, using a small number, 
often just one, of different types of objects. When we are covering the floor of a room 
using identical rectangular tiles (in the ordinary sense of the word) we lay the copies of 
the tiles in a regular way next to each other, leaving no gaps. 

Your floor can be as boring as the rectangular floor on the left or as interesting as the 
Escher lizard tiling on the right: 




In both these floor tilings there is only one tile used, a square tile on the left and a lizard 
shape tile 
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on the right. 

Beyond aesthetics and decoration of the tile there is another difference between these 
two examples which is very significant from our point of view. To fill the floor using the 
square tile we only need to translate the tile around; we never have to turn it. This is not 
the case with the lizard tile. To fill the floor the various copies of the lizard tiles have to 
be turned appropriately so as to fit perfectly with each other. 

It is with the seemingly boring case that we shall occupy ourselves in this survey: the 
single tile that we have at our disposal is only allowed to translate in space. We cannot 
rotate it or reflect it. True, in the vast majority of examples in the literature (see e.g. 
[13]) that exhibit interesting behavior (such as undecidability or aperiodicity, on which 
we'll have more to say later) the allowed group of motions is usually the full group of rigid 
motions in space and therefore rotations are allowed (and some times even reflections as 
well) . One would indeed be hard pressed to find eye-catching examples such as the Penrose 
"kite and dart" aperiodic set of tiles shown below. 




But we have good reasons for restricting ourselves to tiling by translation which we hope 
to expose in this survey. To begin with let us state that tiling by translation is still full of 
unresolved questions with connections to number theory, Fourier analysis and the theory 
of computation, and that happens even in dimension 1. If it is hard to imagine how one 
can seriously study questions of tiling in dimension 1 let us point out that the tiles need 
not have a "nice" shape. For instance, one can easily show by trial and error that the set 
E = {0, 2, 3, 5} (show below in red) 

• o • • o • oooooo 

can tile the integers Z by translation only and that such a tiling has period 12. This trial 
and error process can become slower and slower if one considers larger finite sets of integers 
and simply asks if they could be used to tile Z. In fact, it is not known how to do this in 
essentialy any other way. This reflects how little we understand this tiling phenomenon. 
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Here is another example, of a more geometric flavor, to convince the reader that taking 
away the freedom to rotate and reflect the tile does not take away the fun. For this consider 
the so called notched cube (left) 




and try to visualize a way to fill 3 dimensional space with it (here it would not help to rotate 
or reflect it). We are making no assumptions here about the lengths of the rectangular cut 
(notch) at one corner of the cube or rectangle. Whatever these lengths are this notched 
cube can indeed tile space by translation [18]. Yet the simplest way we know how to prove 
this is using Fourier analysis, making absolutely no use of geometric intuition. Same holds 
for the extended cube (above right) and this seems to be even harder to visualize. At least 
the notched cube "forces" a corner of a copy of the tile to fill the notch in another copy 
and this can at least get you started but no such restriction is obvious for the extended 
cube. 

By now we should fix the rules of the game. What does it mean for a subset Q of W 1 to 
tile when translated at the locations A C 18L d (another discrete set)? Simply that the 
copies 

Q + A, A 6 A, 

are mutually non- overlapping, and cover the whole space. Here, as is customary, we denote 
by O + A the set 

{lu + A : we^}, 




or, in other words, the translate of f2 by the vector A. And what does it mean for two 
sets to be non-overlapping? One could, for instance, demand that the interiors, but not 
the boundaries, of the two sets are disjoint. But it turns out that the most convenient way 
is to demand that the intersection of these two sets has zero volume (or area in dimension 
d = 2, or length in dimension d = 1) or, more precisely, zero Lebesgue measure. 

2. Using the Fourier Transform 

2.1. Tiling in the Fourier domain. Having specified what we mean we turn now to 
exploiting the technology that we know, and in this case we find it very fruitful to define 
tiling simply by the equation 

(1) Xfi( x ~ A) = 1, (almost everywhere for x € 

xeA 
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We now have an equation! as a physicist would exclaim. A mathematician should never 
underestimate the power of formal manipulation and the phenomena it could reveal. So, 
let us rewrite our equation (ignoring from now the exception of measure): 

(2) xn * <5a = 1- 

Here the * operator denotes convolution and 5\ is the measure 




<5a = ^ 8\, 

AsA 

where 5\ is a unit point mass sitting at point A € M. d . (We recall that / * 5\(x) — f(x — A), 
so that convolving a function / with 5\ merely translates the function by the vector A.) 
The object 5\ is a compact way to encode all information that's contained in the set of 
translates A in a way that we can operate algebraically on it. 

An equation such as (2) of course begs to be subjected to the Fourier Transform 1 as 
the Fourier Transform (FT) behaves so nicely with convolution (the FT of a convolution 
is the pointwise product of the FTs of the convolution factors: fg = fg). The definition 
of the FT / of a function / that we use is 

(3) /(£)=/ f(x)e- 2 ^*dx, (eG< / |/| < oo). 

Jm. d Jm. d 

This definition is sufficient for the function / = xn which is integrable but one needs to 
define the FT of a measure by duality (see for instance [22] for the details). Postponing 
such worries till later we take the FT of both sides of (2) to get 

(4) • 5 A = S . 

The measure So on the right is just a point mass at point 0. The support (where it is 
"nonzero") of this is {0} and therefore that must be the support of the left hand side of 

(4) . Since that is a product it follows that, apart from point 0, wherever xn is nonzero, S\ 
must be zero so as to kill the product. Let us summarize this in the following necessary 
(and often sufficient) condition for tiling: 

(5) supp^C {xh = 0}U{0}. 

A formal proof of (5) can be looked up in [22]. Let us make no comment here about when 
this condition is sufficient. But we must point out that the object S\ is not necessarily a 
measure but a so called tempered distribution [35]. 

2.2. A problem from the Scottish Book. Nice reformulation but what could this buy 
for us? 

Following the time-honored tradition of our trade, let us first generalize: whatever we've 
said so far holds just as well for an (almost) arbitrary function / (nonnegative, integrable 
would do) in place of xn- In other words, if one has 

(6) ^ f(x - X) = £ = Const., 

AeA 

for almost every x 6 K d then it follows that 

(7) supp^C Z(/)U{0} where Z(/) = {/ = o}. 



This of course depends on what training you have been subjected to. 
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Whenever (6) holds we will say from now on that / tiles space with A at level £. And 
Steinhaus suggested that perhaps f(x) = e~ x has this property (although he quickly 
disproved it). 

Let us apply our new tool to Problem No 181 in the famous Scottish Book 2 [31], a 
problem posed by H. Steinhaus and partly solved by him around 1939. 



(in the sense of H. Steinhaus) for every couple f,,/,(f, * U)l 


181 


Find a continuous function {or perhaps an analytic one) 
f{x), positive and such that one has 

OP 


I V ■ 

H. STEINHAUS 


E f{X + B) ■ 1 
ft - - » 




(identically in x in the interval -a> < x < + »); examine 
whether {\/-^e~'' is such a function; or else prove the 
impossibility; or else prove uniqueness. 




Addendum. The function (l/VVJe"*' does not have the 
property — this follows from the sign of the second derivative 
For v = of the expression 




E -U-'-"' 1 . 

- v H. Steinhaus 





In our language the question essentially was if there is a positive, continuous and even 
analytic function / that tiles M with Z. 

2 

Whoever has been exposed even a little to the FT knows that e~ x has the nice property 
that its FT transform is essentialy itself (times a constant and possibly rescaled, depending 
on the normalizations used in the definition of the FT) and therefore it has no zeros at 
all and therefore (7) has no chance of holding as <5a has no place to be supported at apart 
from 0, and it's not a constant. And this argument does not use the fact that A = Z so 

2 

e~ x does not tile M with any set of translates. 

An analytic solution to this question of Steinhaus can also be provided using our tech- 
nology. First of all, let us make the important remark that 

(8) 6 Z d = 6 X d . 

This is the deservedly famous Poisson Summation Formula [36, Ch. 9] 

n n 

in disguise. In this case, when 6a is a measure, condition (7) is also sufficient for tiling [22], 
so it is sufficient to find a positive function / whose FT vanishes on Z \ {0}. Analyticity 
of / will be guaranteed if / has compact support. It is then easy to check that one can 
take as / the sum of two triangle functions with incommensurable base lengths, 



as the FT of each such triangle is a nonnegative function (a Fejer kernel) whose zeros are 
at integer multiples of the reciprocal of its half-base length. Since these will never match 

2 A collection of problems from the Scottish Cafe in Lwow (then in Poland now Lviv in the Ukraine) 
where a number of important mathematicians met for a number of years exchanging problems. These 
problems and whatever solutions were recorded in a notebook that formed the basis for [31]. 
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for the two Fejer kernels their sum is always positive. Since we also take the bases of the 
two triangles to be supported in ( — 1, 1) it follows that (7) holds and we have tiling. 

2.3. Filling a box with two types of bricks. Let us give another amusing application 
of the FT method to tiling problems. This will be a Fourier analytic proof [21] of a result 
of Bower and Michael [5] . Suppose that you have two types of rectangular bricks at your 
disposal, type A of dimensions a\ x (12 and type B of dimensions b\ x 62 (we're stating 
everything in dimension 2 but everything works in any dimension) and your task is to tile 
a rectangular box, say Q — (—1/2, 1/2) 2 , with copies of bricks A and B. The bricks may 
be translated but not rotated. 

We will show that this is possible if and only if you can cut Q along the x or along the 
y direction into two rectangles each of which can be tiled using bricks of one type only. 
A generic tiling of the type shown below left implies the existence of a separated tiling as 
shown below right. 





In other words, having two types of bricks at your disposal does not demand much ingenuity 
on your side in order to exploit for tiling purposes. 

A simple calculation shows that ifC = ( — y, y) x (— %) is a centered c\ x C2 box 



then 



3in(7rci£) sin(7TC2?7) 



and therefore xc vanishes at those points, shown below, 











V 
































l/c 2 




































1/d 







































where the £ or the r/ coordinate is a non-zero multiple of 1/ci or I/C2 respectively 

Suppose now that we can tile Q by translating copies of brick A to the locations T and 
copies of brick B to the locations S. In other words 



(9) 

which we rewrite as 



XQ = XA* St + Xb * Sg 
and take the FT of both sides to get 

(10) xq(£, v) = <h(€, v)xa(£, v) + 0s(£. v)xb(^ v), 
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where (fix = <5r and <f>S — °S are two trigonometric polynomials. Since Q — (—1/2, 1/2) 2 
we have 

(11) z(xq) = {x^ = o} = {Zez\{o}or V ez\ {o}}. 

But, because of (10), \Q must vanish on the common zeros of \X and xb, for instance at 
the points (1/a.i, I/&2) and (l/&i, 1/0-2)) which implies, because of (11) 

(12) (l/oi 6 Z or l/6 2 6 Z) and (l/&i G Z or l/o 2 G Z). 

If (12) is true because l/oi, l/a 2 6 Z then brick A alone can fill Q. Similarly if (12) is 
satisfied with l/&i, I/&2 6 Z then brick B alone suffices. 

What happens if l/oi, 6 Z? Since we have assumed a tiling of Q in (9) it follows, 
by traversing the y-axis, that there are nonnegative integers k, I such that 

1 = ka 2 + lb 2 . 

Cut now the Q box parallel to the x-axis at height ka 2 from the bottom as show here: 



1 













ka 2 



Now it is clear that brick A can fill the lower box (since 1 ja\ G Z) and brick B can fill the 
upper box (since l/&i G Z). The remaining case l/a.2, I/&2 G Z is treated similarly. 

3. Discrete tilings 

3.1. Tilings of the integers and periodicity. Let us now focus on tiling the integers. 
Let A C Z be a finite set and A G Z. We say that A tiles Z with A at level I if the copies 
A + A, A G A cover every integer exactly I times. In other words 

xaOe — A) — £, for all x G Z. 

AeA 

We denote this situation as A + A — PL. 

We say that a tiling is periodic with period t G Z if A + 1 — A. It is a basic fact proved 
by Newman [32] that all tilings of the integers at level 1 are periodic. 

Indeed, suppose that A = {0 — a\ < . . . < a^} is a finite set (we may freely translate A 
without changing the tiling or the periodicity property so we assume it starts at 0) and 
A + A — Z is a tiling at level 1. Fix any x G Z and write W x = x, x + 1, . . . , x + — 1 
for the "window" of width afc — 1 (one less than A) starting at x. We claim that the set 
A fl W x determines A. Let us show that it determines A to the right of x + — 1. It will 
determine A to the left of x by the same argument. 

It is enough to decide, looking at A n W only, if x + cik G A or not. (We then repeat for 
x + dk + 1 and so on.) Observe that for any A G A fl (—00, x) the set A + A is contained 
in (—00, x + dk — 1], so any such copy cannot be used to cover x + a*;. Clearly no copy of 
the form A + A for A > x + at can be used for that purpose too. We conclude that x + 
is covered by some copy A + A with A G W U {x + a^}. Inspecting A fl W x we can tell if 
the relevant A is in W x or not. If it is in W x then x + £ A since this would lead to the 
copies A + A and A + (x + at) to overlap at x + a^. If it is not in W x then necessarily 
x + afc G A, and this concludes the proof of the claim. 

How many different values can the set A fl W x take? Clearly it can take at most 2 afc 
different values as there are two choices (in A or not in A) for each x G W x . This means 
that there are two different x, y G {0, 1, . . . , 2 afc } for which A(~)W X is a translate of AnW y . 
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It follows that A + (y — x) = A. We have proved that every tiling has a period which is at 
most 2 D where D is the diameter of the tile. 

There is a similar result for tilings of the continuous line by translates of a function 
[23, 28, 22]. Combinatorial arguments do not seem to be enough here and the Fourier 
analytic technology along with some deep results of Harmonic Analysis are being used for 
the proof. 



3.2. Tilings of the finite cyclic group. The fact that a tiling of the integers is periodic 
allows us to view it as a tiling on a smaller structure, a cyclic group. Indeed, assume that 
the tiling A + A = Z has period n, that is A + n — A. Define the set 

A — A mod n C {0, 1, . . . , n — 1} 

by taking for each A £ A its residue mod n. It follows from the n-periodicity of A that 
A — A + nZ (this is again a tiling or a direct sum: every element of A can be written in a 
unique way as an element of A plus an element of nZ) . Hence we have 

Z = A + K = A + K + nL 

with all sums being direct. Taking quotients we obtain that the cyclic group Z„ — Z/(nZ) 
can be written as a direct sum (tiling) 

(13) Z„ = A + A. 



In this case we obviously have n — \A\ 



A 



Let us stop here to make two side remarks about periodicity. The first is that tilings 
of the cyclic group Z„ can also be periodic. Indeed, assume Z n — A + A and there exists 
/ k £ Z„ such that A is periodic by k, i.e. A + k — A. Then we can reduce the set 
A modulo k, A = A mod k C {0, 1, . . . , k — 1}, and conclude as above that Z& — A + A. 
Therefore, periodic tilings of Z n can be regarded as tilings of a smaller group (it is 
trivial to see that k automatically divides n). It is thus natural to ask whether certain 
cyclic groups Z n admit only periodic tilings, i.e. whenever Z n = A + A, then either A or A 
must be periodic. These were called "good" groups by Hajos [14] (the notion also makes 
sense in the more general setting of finite Abelian groups, not only cyclic groups). It turns 
out that some groups indeed have this property, and Sands completed the classification 
of good groups in [38, 39]. In particular, the good groups that are cyclic are Z„ where n 
divides one of pqrs, p 2 qr, p 2 q 2 or p m q, where p, q, r, s are any distinct primes. The cyclic 
group of the smallest order which is not good is Z72. 

The other remark concerns the connection of tilings of cyclic groups to music com- 
position. As explained above, periodic tilings of Z„ are mathematically speaking less 
interesting because they can be considered as tilings of some smaller group Z&. It turns 
out that non-periodic tilings of a cyclic group Z„ are also more interesting from an aes- 
thetic point of view, and they are called Vuza- canons in the musical community. The 
interaction of mathematical theory and musical background has been extensively studied 
in recent years [3, 1, 2, 10, 42]. Finding all non-periodic tilings of Z„ is therefore motivated 
by contemporary music compositions. Of course, the problem makes sense only if Z™ is 
not a good group (otherwise all tilings are periodic). Fripertinger [9] achieved this task for 
n = 72, 108 while the authors [26] gave an efficient algorithm to settle the case n — 144 
(the algorithm is likely to work for other values like n — 120, 180, 200, 216, beyond which 
the task simply seems hopeless.) 

Let us now return to the analysis of tilings of the cyclic group Z„. The Fourier condition 
(5) for tiling takes exactly the same form here and is much simpler to prove as no subtle 
analysis is required (no integrals, only finite sums are involved). The Fourier transform of 
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a function / : Z n — > C is denned as the function / : Z n — > C given by 

n-l 
3=0 

where £„ — e 2m / n is a primitive n-th root of unity. 

In the finite case the roles of the tile A and the set of translations A in (13) can be 
interchanged, so let us adopt a more symmetric notation 

Z n = A + B, 

in which A, B C Z„ are just two subsets, necessarily satisfying \ A\ ■ \B\ — n, and such that 
every element of Z„ can be written uniquely as a sum of an element of A and en element 
of B. The Fourier condition now takes the form 

(14) Z n = {0}UZ(x2)UZ(x^). 

As we saw in the section §3.1 every tiling of Z has a period which is at most 2 D , where D 
is the diameter of the tile. We will see in this section how the easy combinatorial argument 
of §3.1 can be replaced with an argument that has the Fourier condition (14) as a starting 
point, uses some well-known number-theoretic facts (cyclotomic polynomials) and gives 
much better results. We will describe the result in [20] and by I. Ruzsa in an appendix in 
[41] since it is simpler than the current best results in [4]. 

So, suppose that A C {0, 1, ... , D} is a set of integers of diameter < D and that A tiles 
Z with period M. This implies that A (to be precise, A reduced mod M) tiles the cyclic 
group Z M 

Z M = A + B. 

Assume also that M is the least period. This implies that if g £ %u is such that B = B + g 
then g = 0. (Otherwise the original tiling would be periodic with period g.) 

Now we use the Fourier transform without really mentioning it, using instead polynomial 
terminology. Let 

A(x) = J2* a , B(x) = Y,x\ 

aeA beB 

be the polynomials defined by A and B. It is very easy to see that the tiling condition 
Z M — A + B is the same as 

A{x)B{x) = 1 + x + x 2 + • • • x M_1 mod x M - 1, 

or, in other words, that 

r M _ i 

(15) x M - 1 divides A(x)B(x) , 

x — 1 

and this implies that all M-th roots of unity except 1 are roots of the product A{x)B(x). 
(This last statement is the equivalent of (14) in algebraic language.) 

The M-th roots of unity ( 3 M = e 2m i/M ^ = 0, 1, . . . , M — 1, are grouped into cyclotomic 
classes. Two roots C,° M , C,\j belong to the same cyclotomic class if and only if the greatest 
common divisors (j, M) and (k, M) are the same. If a root of unity in a certain cyclotomic 
class is a root of a polynomial with integer or rational coefficients then so are all the other 
roots in the same cyclotomic class. If d is a positive integer then the <i-th cyclotomic 
polynomial is the monic polynomial which has as roots all the primitive d-th. roots of unity 

*d[x)= J] (x-e 2 ^ d ). 
l<j<d,(j,d)=l 
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We have deg^^x) — (f>(d) (the Euler function) and we will use the estimate [15, p. 267] 

(16) C r^~^ ^ ^ ^ d ' 

log log a 

which holds for some constant C and all large enough d. It turns out that each is an 
inrreducible polynomial with integer coefficients and one can write 

x m - 1 = n 

d divides M 

Coupled with (15) this implies that if d > 1 is a divisor of M then ^(x) divides the 
product A{x)B(x). 

Let now & si (x), . . . , <& Sk (x) be all cyclotomic polynomials <3> s (:r) with s > 1 that divide 
A(x), written once each and numbered so that 1 < si < s 2 < ■ ■ ■ < Sfc. Since deg $ s — <f>(s) 
it follows that 

(17) 0(si) + • • • + 0(s fc ) < deg < D. 
We have 



X! Si - C^2<j){si) log log s t (from (16)) 



i=l i=l 
fe 

<C^</>( Sl )loglogD 

i=l 

<CL>loglogD (from (17)). 

From this inequality and using the fact that the Sj are different integers it follows that 
(18) k < Cy?D log log D, 

as obviously the worst case is if the Sj are k consecutive integers and in that case their 
sum is > k 2 . 

Define now t = Yli=i s «' so tnat au cyclotomic polynomials that divide A(x) are also 
divisors of x l — 1. It follows that 



,M 



x 



1 divides {x l - l)B(x) 



and this means precisely that B has period t, hence t > M as we assumed M to be 
minimal. Using (18) and the bound Sj = 0(D 2 ), for instance, we get the following bound 
for the least period M: 

(19) M <t< exp(C/DlogD^loglogL»), 

for an appropriate constant C. This is a much better bound, in terms of its dependence 
on D, that the bound M < 2 that we obtained with combinatorics alone. 

3.3. Periodicity in two dimensions and computability. In contrast to the one- 
dimensional case, where every tiling by a finite subset of Z is periodic, it is very easy 
to see that in Z 2 there are tilings which are not periodic. One has to be a little careful 
with periodicity in dimension 2 and higher though. We shall call a set A C 7L d periodic 
if there exists a full lattice of periods, i.e. if there exist d linearly independent vectors 
ui,...,Ud which are all periods of the set 

A + u\ — A + u 2 = ■ ■ ■ = A + u d — A. 

If that happens then every integer linear combination of the Uj is also a period or, in other 
words, the set A has a full-dimensional lattice of periods 

A = span z {u l ,u 2 ,...,u d }. 
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Thus in dimension 2 a set may have a period but not two linearly independent periods, in 
which case we do not call it periodic. One example is the set Z x {0} C Z 2 . A tiling of Z 2 
by a finite subset which is not periodic is very easy to construct. 



On takes A to be a square, for instance the set {(0,0), (1,0), (1,0), (1, 1)} and perturbs 
the usual tiling of it by shifting one column only up by one unit. This destroys the period 
(2, 0) of the tiling along the x-axis, but leaves the period (0, 2) along the j/-axis intact. 

It is not entirely obvious how to construct a tiling of Z 2 which has no period at all, 
but it can be done and let us briefly describe how. The key idea is to start again with 
the usual square tiling and do something to it so as to destroy all periods, not just the 
periods along one axis. The way to do this is to simultanesouly shift a horizontal and 
a vertical column. This is not of course possible with the square tile we used before, as 
simple experimentation will convince you. But it can be done if we "interleave" four tilings 
of the type shown above. Take our tile to be the set 

A = {(0,0), (2,0), (0,2), (2,2)} 

which is what we had before only dilated by a factor of two. Now tile (2Z) 2 using A in 
the usual way, that is by translating A to the locations (4Z) 2 . Using this tiling one can 
get a tiling of Z 2 by tiling the four sets (cosets of the subgroup (2Z) 2 in Z 2 ) 

(2Z) 2 , (2Z) 2 + (1,0), (2Z) 2 + (0,1), (2Z) 2 + (1,1) 

in exactly the same way. But the smart way to do it is to use the nice tiling for the first two 
of the above four cosets, then use, for tiling (2Z) 2 + (0, 1), a tiling like that shown in the 
figure above (which destroys the horizontal periods) and finally use for tiling (2Z) 2 + (1, 1) 
a similar tiling which destroys the vertical periods. 

Having established that in two dimensions there are tilings with no periods at all, let us 
now remark that it is very different to ask for tiles which are aperiodic, i.e. for tiles which 
can tile but only in a way that admits no periods. In fact, the answer to this question is 
not known if we insist that only translations are allowed. 

Conjecture 3.1 (Lagarias and Wang [28]). If a finite subset A C Z 2 can tile by translation 
then it can also tile in a periodic way. 

Let us point out that if more freedom than translation is allowed then tiles (or sets of 
tiles) are known which can only tile aperiodically, so the above conjecture, if true, would 
mean that restricting the allowed motions to translations makes a big difference in this 
respect. 

Let us also remark that it can be proved, at least in dimension 2 [33], that if a set 
admits a tiling with just one period then it also admits a (fully) periodic tiling. 

The property of periodicity is of great importance to the issue of computability, namely 
to whether a computer can decide if a finite set A admits tiling by translation. In still 
more words, we are interested to know if there is a computer algorithm (a Turing machine 
for purists) which, given a finite subset A oil? 1 will decide in finite time if there exists a 
tiling complement of A or not, i.e. if there exists B C Z 2 such that A + B = Z 2 is a tiling. 

We do not require our algorithm to: 

• run fast (only to finish at some point) 
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• find such a tiling complement B if it exists. 

Especially for the second point notice that it would not make sense to want to know B as 
such a set is an infinite object and there is no a priori reason for it to be describable in 
some finite manner. 

It is easy to give an algorithm that would answer NO if A is not a tile but would 
run forever if A is a tile. Here is how this can be done. For each n — 1,2,..., decide 
by trial and error if there is a finite set B such that A + B covers the (discrete) square 
Q n = [— n,n] 2 n Z 2 in a non-overlapping way. There is a finite number of such sets B to 
try as it does not matter what B is far away from Q n . If such a B exists move on to 
the next n. If not then declare A a non-tile and stop. The correctness of this algorithm 
(i.e., the fact that the algorithm will stop for any non-tile A) follows from a fairly simple 
diagonal argument: if A can cover in a nonoverlapping way every Q n then it can also tile 
Z 2 . We invite the reader to think over why this is so. 

The hard thing of course is to have an algorithm that always halts and say YES if A is 
a tile and NO if it is not a tile. Let us point out here that in various other tiling situations 
such an algorithm cannot exist. For instance [12] there cannot exist an algorithm which, 
given a finite collection of subsets of Z 2 , decides if these can tile the plane by translation. 
This remains true even if the number of subsets remains bounded (but large) . 

The connection with periodicity [34] is that the periodic tiling conjecture above implies 
decidability. Indeed, let us assume that any set that tiles can also tile periodically. (We 
emphasize here that we are making no assumption about the size of the periods or their 
dependence on the tile.) This means that it is equivalent to decide if a given set A has 
periodic tilings. 

Let us assume that A C [0, D] 2 . For n — 1,2,..., we can clearly find all subsets B of 
[—n — 2D, n + 2D] 2 such that A + B is a nonoverlapping covering of Q n — [— n, n] 2 . We can 
do this, slowly but surely, by just trying all eligible subsets. If we find none then clearly 
our set is not a tile, we say NO, and we stop. If we do find some we look among them to 
find a set B which "can be extended periodically" . How does such a set look like? 

If an infinite set B C Z 2 is periodic with linearly independent period vectors U\,U2 6 Z 2 
then we can always find two vectors 

ui = (a, 0), u 2 = (0, b), a, b £ {1, 2, 3, . . .}, 

which are also periods (in other words, every lattice in Z 2 contains a sublattice generated 
by two vectors on the x- and y-axes). The mapping (x, y) — > (x mod a, y mod b) maps B 
to a set 

B' C {0, 1, . . . , a - 1} x {0, 1, . . . , b - 1} 
and the periodicity of B means exactly that 

B = (aZ) x (bZ) + B' 

is a direct sum, i.e. every element of B can be written uniquely in the form (am, bn) + b' , 
where m, n 6 Z, b' 6 B' . 

Suppose now that our algorithm looks at a finite set B which is a finite part of a periodic 
set B. This finite part of B arises if we keep those elements of B which are used in a non- 
overlapping covering of Q n (shown as a dashed green rectangle). If n is large enough then 
the set B will look something like this: 
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That is, the set B will contain several "full periods" (the dashed blue a x b rectangles) 
B plus a few incomplete periods near the border of Q n . Suppose now that we have a k x k 
block of full periods in the set B and that k is large enough that both ka and kb are much 
larger than the diameter of A, say larger than 10D. Then we can just extend this block 
periodically (with periods (a, 0) and (0, b)) and we will get a periodic tiling complement 
of A. The reason is that (a) no overlaps arise in this manner or they would have shown 
up in the k x k block already and that (b) we have covering of everything as clearly the 
k x k block suffices to cover the square of side 5D with the same center. 



4. Tiles, spectral sets and complex Hadamard matrices 

In this last section we introduce the notion of spectral sets and describe how it is related 
to tilings and complex Hadamard matrices. For the sake of simplicity, we will restrict our 
attention to finite groups, in particular to cyclic groups Z ra and its powers Z^. 

Although we have seen in the sections above that many fundamental problems about 
translational tilings remain open, the reader will agree with us that the notion of tiling 
is very intuitive and easy to grasp. To the contrary, the definition of spectral sets is 
somewhat more abstract. In order to make it more down-to-earth, we will not use the 
standard definition here, but first introduce complex Hadamard matrices and then use 
them to define spectral sets. 

What are Hadamard matrices? Classically, they are square matrices consisting of el- 
ements ±1 only, such that the rows (and thus also the columns) are orthogonal to each 
other. Complex Hadamard matrices are a natural generalization of this concept. A k x k 
matrix H is a complex Hadamard matrix if its entries are complex numbers of modulus 
1, and the rows (and thus the columns) are orthogonal. Recall here that orthogonality is 
understood with respect to the complex scalar product, i.e. you need to conjugate in one 
of the components, the scalar product of (z\, . . . Zk) and (m, . . . Uk) being Ylj=i ~Zj u j- 

It is well-known that a k x k real Hadamard matrix can only exist if k is divisible by 4. 
On the contrary, complex Hadamard matrices exist in all dimension. Indeed the Fourier 
matrix Fk defined as 

11 1 

1 Ck 



F k := 



/■k-1 
^k 
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where (jfc — e 2m / k , is a complex Hadamard matrix for every k. 

Each element hj_ m of a k x k complex Hadamard matrix H is of the form e 2 ™ p ^ m where 
Pj,m € [0, 1)- We will call the real k x k matrix R formed by the angles p^ m the logarithm 
of H, in notation, R = \og{H) or H = exp(i?). For a real matrix R such that exp(i?) is 
complex Hadamard, we will use the terminology that R is a log- Hadamard matrix. 

Complex Hadamard matrices play an important role in quantum information theory, 
in particular in the construction of teleportation and dense coding schemes [43] , among 
other things. An online catalogue of all known families of complex Hadamard matrices is 
available at [6]. 

Let us now turn to the definition of spectral sets. We identify elements of Z^ with column 
vectors, each coordinate being in the range {0, 1, . . . , n — 1}. Thus, a set S C with r 
elements can be identified with andxr matrix, the columns of which are the elements of 
S (the order of elements does not matter) . We will abuse notation and denote this matrix 

also by S. We also identify Z^ with row vectors whose coordinates are in {0, 1, . . . , n — 1}. 
Accordingly, we will use the notation Z^ — Z^ T (the T meaning transposition). Two sets 
S C Z^ and Q C Z^ T , both of them having r elements, are called a spectral pair if the 
matrix ^QS is an r x r log-Hadamard matrix. In this case Q is called the spectrum of S, 
and S is called a spectral set. We remark that S T Q T = (QS) T (where S T and Q T denote 
the transposed matrices) is then automatically also a log-Hadamard matrix, so we obtain 
that S is a spectrum of Q, which justifies the symmetry of the terminology "spectral pair" . 

For example, the whole group Z n = {0, 1, . . . , n — 1} is spectral, and its spectrum is the 
whole dual group Z^, giving rise to the logarithm of the Fourier matrix, log(F n ) = -Z^Z n . 

We saw that the tiling Z„ — A + B is characterized by \A\ ■ \B\ = n and the Fourier 
condition Z n = {0} U Z (xa) U Z (xb) in equation (14). After some inspection one sees that 
spectral sets also admit a characterization by Fourier analysis. Indeed, using the definition 
above we obtain that S C Z^, \S\ = r, is spectral if and only if there exists Q C Z^ T , such 
that |Q| = r and 

(20) qj-q k eZ(x~s) 

for each qj / 6 Q- (This is an equivalent way of saying that the jth and kth rows of 
the matrix exp(iQS') are orthogonal.) 

What are the connections of tiles and spectral sets? Notice that by the orthogonality 
conditions the matrix U — ejcp(^QS) is unitary. This implies the following important 

fact (for a formal proof, see [25]): 

(21) ^|xS| 2 (*- g ) = |S| 2 = r 2 , 

q&Q 

for all x G Z^ T . In other words, the function \xs\ 2 tiles Z^ at level \S\ 2 when translated 

at the locations Q. One can also say that the set Q tiles the group Z^ with the weighted 
translates defined by the nonnegative function 




The question thus arises naturally: does Q tile Z^ in the ordinary sense, i.e. does there 

exist a set P c Z^ such that Q + P = Z^? A famous conjecture of Fuglede [11] concerns 
exactly this scenario: it states that a set S is spectral if and only if it tiles the group 
(the conjecture was originally stated in the Euclidean space R d but it makes sense in any 
abelian group). 

Several positive partial results were proven concerning special cases of Fuglede's conjec- 
ture (see e.g. [16, 17, 19, 27]), before T. Tao [40] showed an example of a spectral set in 
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M 5 which did not tile the space. Tao's example was based on considerations in the finite 
group Zjj. It also implies the existence of such sets in any dimension d > 5. Subsequently, 
counterexamples in lower dimensions (4 and 3, respectively) were found by the authors 
[29, 24]. All these examples are based on considerations in finite groups Z^ and, ultimately, 
on the existence of certain complex Hadamard matrices. 

The other direction of Fuglede's conjecture could not be settled by Tao's arguments. 
Finally a non-spectral tile in IR 5 was exhibited by the authors [25] by a tricky duality argu- 
ment, also based on considerations in finite groups. Later, counterexamples in dimensions 
4 and 3, respectively, were found in [8, 7]. As of today, Fuglede's conjecture is still open 
in both directions in dimensions 1 and 2. 

As a final remark let us mention another interesting connection of tilings and complex 
Hadamard matrices. As discussed above, Fuglede's conjecture is not true in general, but it 
is true in many particular cases. What could this buy for us? We saw from the definition 
that spectral sets are directly related to complex Hadamard matrices, and the latter are 
very useful in another branch of mathematics: quantum information theory. Could we use 
the connection of tilings to spectral sets (which exists in many special cases) to construct 
new families of complex Hadamard matrices? It turns out that the answer is positive. It 
turns out [30] that "trivial" tiling constructions unfortunately lead to well-known families 
of complex Hadamard matrices (so-called Dita- families) . However, a non-standard tiling 
construction of Szabo [37] was used in [30] to produce previously unknown families of 
complex Hadamard matrices in dimensions 8, 12 and 16. 
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